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Conventional many- body calculations which do not consider mesonic degrees of free- 
dom are not reliable for the study of nuclear matter where high densities are present. For 
these purposes, a relativistic calculation must be performed. Quantum chromodynamics 
(QCD) is the fundamental theory of the strong interaction and hence it should explain 
possible modifications of hadron properties in the nuclear medium. However, typical nu- 
clear phenomena at intermediate and low energies cannot be analitically derived from 
QCD although one hopes that QCD will be solved numerically on the lattice in a near 
future. Meanwhile we are left with the construction of phenomenological models in order 
to try to describe nuclear phenomena and bulk properties. Walecka and collaborators 
used such a kind of theory for the first time around 1974 to describe the nucleon - nucleon 
interaction [Q. Since then the Walecka model has been widely used to describe the 
properties of nuclear matter as well as some properties of finite nuclei. This model is 
based on a phenomenological treatment of the hadronic degrees of freedom and, because 
of this fact, it is also known as QHD-I (quantum hadrodynamics) , which consists in a 
renormalizable relativistic quantum field theory. The early version of this model consid- 
ers a scalar (ex) meson field and a vector (o>) meson field coupled to the baryonic (nucleon) 
field. Besides the relativistic mean field calculation, the Walecka model has also been used 
in a more complete treatment, the relativistic Hartree- Fock approximation ||, ||, ||. 

Some of the drawbacks of the model are that the effective nucleon mass obtained at 
high densities is too small and its incompressibility at the energy density saturation is 
too large. To avoid this problem, Zimanyi and Moszkowski [||] introduced an alternative 
coupling between the scalar meson and the nucleon (which is a scalar derivative coupling) 
and another coupling between the scalar and vector mesons. We should stress that a 
microscopic foundation for derivative scalar coupling models has been derived from the 
relativistic SU(6) model |J. It has already been investigated, in a mean- field approx- 
imation @,||, how this derivative coupling Zimanyi - Moszkowski (ZM) model differs 
from the usual Walecka model when the effective nucleon mass, the energy density, the 
incompressibility and other important features are calculated. 

To try to solve the above mentioned drawbacks of the Walecka model, more sophisti- 
cated treatments have also been developed. One of them || refers to a relativistic Hartree 
- Fock approximation in a non - linear Walecka model for nuclear matter and finite nuclei 
which considers a, u, n and p mesons. In this work, the a - meson mass is density depen- 
dent and the non linear character of its related Hamiltonian makes the calculation of the 
exchange contributions to the energy rather complicated, which requires an approximation 
for the calculation of this Fock - like terms. Even though, the inclusion of the exchange 
terms and the isovector mesons (n, p) turn out to be important for the description of 
the energy and the incompressibility in nuclear matter and of the spin -orbit interaction. 
Another successful treatment utilizes the so called Relativistic Dirac-Brueckner-Hartree- 
Fock method (RDBHF ) [|T^]. In this case the coupling constants for the sigma and omega 
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mesons become density dependent and they actually decrease with increasing energy. The 
computational effort also increases considerably since self - consistency for the self ener- 
gies is obtained by considering different coupling constants at each density. Nevertheless 
this sort of calculation has proved very good when applied to finite nuclei and the 



authors concluded that Fock terms are small but important and should not be neglected. 

It is important to stress that the ZM models also belong to the category of descrip- 
tions involving density-dependent coupling constants ||. The description underlying the 
approach known as relativistic density- dependent Hartree-Fock [12| reproduces finite nu- 



clei and nuclear matter saturation properties using coupling constants that are fitted, at 
each density value, to the RDBHF self-energy terms. The good agreement obtained for 
the ground state properties of spherical nuclei lends support to this sort of description 
involving density- dependent coupling constants. 

In the present work, we propose an alternative way of considering in a very simple 
manner direct and exchange terms in the infinite nuclear matter to obtain the properties 
of the ZM models mentioned above and compare them with the results which come out 
of the relativistic Hartree - Fock approximation for the Walecka model. Our aim is to 
use these derivative scalar coupling models to obtain bulk properties in finite nuclei in a 
forthcoming work and for this purpose, based on the conclusions drawn in refs. [H, |§, 



11 1, we assume that the exchange terms will be necessary. In this sense, the present work 
has also a didatic meaning. 

In the recent literature three different possibilities have been considered for the cou- 
pling of the nucleon with the mesons in the Zimanyi - Moszkowski (ZM) model (for a 
review, check ref. 0). In this letter we consider two of them. The first one is known 
as the original ZM model and the second one has been chosen because is the one which 
gives better results for the nuclear matter in the mean field approximation and we call it 
ZM3 (to be consistent with the definition in other published material). To start with, we 
write the Lagrangian density for the ZM and ZM3 models respectively as ||: 

C zm = -ijMt/j + kj- 1 ^^ - 9vli y^ 

+ ^(d^dy - m 2 ^ 2 ) - \F, U F^ + im^V, (1) 



and 



£ Z m3 = -tfjMtfj + (m* m ) V^T/^'V - g^^V^tp 
+ l -{d^ - m 2 ^ 2 ) - ^ V + im^V, (2) 



where 



< m = a + ^r 1 (3) 
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and if), 4> and are field amplitudes for the nucleon, the scalar-isoscalar meson, the 
vector-isoscalar meson, F^ u = d^V v — d v V^ and M is the bare nucleon mass. Notice that 
these derivative coupled Lagrangians are Lorentz invariant, but they are not renormaliz- 
able. The physical meaning of this modified couplings is that the kinetic fermionic term 
describes the motion of a particle with an effective mass M* instead of the bare mass M 
present in the conventional Walecka model. 

In ref. ||, a rescaled Lagrangian is obtained from the above equations, with the 
rescaling of the fermion wave function as ip — > y/m* zm ip in eqs.fll]) and (@) and the rescaling 
of the field as — > Tn* zm V^ in eq.(@) and it reads 

C R = $(i7„3" — M + rnZjgdW 

+ i(9 M ^_m^ 2 ), (4) 

where, for the Walecka model a = 0, (3 = 0, for ZM a = 0, (3 — 1 and for ZM3 a = 2, 
>' 1- 

In a very simple description of the ZM and ZM3 models they can be viewed as models 
with non-linear effective scalar coupling constants g* and g* in such a way that the new 
Lagrangians and related quantities can be obtained simply by substituting the old coupling 
constants in the Walecka model by the new ones ||,i.e.: 

C-zm £\V alecka\9 's * 9s) 

and 

where g* = m* m ■ g s and g* = m* m ■ g v . Instead of the above substitutions, which can be 
cumbersome to be carried out in the self energy expressions, we utilize a new subscription 
which introduces modifications in the meson masses only WM : 

C-zm = C-Waleckai^s ~^ m * s ) (5) 

and 

^zm3 = £walecka{m s -> m* s ,m v — ► 771*), (6) 

where m* = m s /m* w and m* = m„/m*y, where, in this case, 



In the Walecka model the relativistic Hartree-Fock equations |H],H are obtained by 
using Dyson's equation to sum to all orders the self-consistent tadpole and exchange 
contributions to the baryon propagator 

G(k) = G°(k) + G°(k)Z(k)G(k), (8) 

where £ is the proper self-energy. Because of translational and rotational invariances in 
the rest frame of infinite nuclear matter and the assumed invariance under parity and 



time reversal, the self-energy may be written as [pi] , [ 13 



E(Jfe) = Z s (k) -7 S°(A;) + ^-kE v (k). (9) 

Coupled integral equations are then solved for the self-energies in the so-called Dirac- 
Hartree-Fock approximation [El. In this approximation just the contributions from real 
nucleons in the Fermi sea are kept in the baryon propagators. The effects of virtual 
nucleons and anti-nucleons on the medium are neglected as usually done in the literature 
0,0. The nucleon propagator in a Fermi sea with Fermi momentum kp is then written 
as (the nuclear density is pb = 2k F /2,n 2 ) 



G(k) = ( 7 ,P* + M *( k ))^ S ( k - E(k)W(k F - |*|), (10) 
k"* = k" + W{k) = (k° + E (Jfe), fc(l + £*(*;))) , (11) 



where 



E*(k) = \J(k*) 2 + M*{kf (12) 
M*(k) =M + E s {k), (13) 
and E(k) is the single-particle energy, which is the solution of the transcendental equation 

E(k) = [E*(k) -E°(k)] h o =m . (14) 

Performing the q° and angular integrals in the expressions for the various components E s , 
E°, S" of the self-energy, three coupled nonlinear integral equations are then obtained. 
For the Walecka model, the equations for the self - energies are given in ref. ||, (pg. 131) 
and we do not reproduce them here. In the ZM model the very same steps are performed 
in order to obtain these three components of the self energy, with the exception that we 
have included an approximation. This approximation amounts to considering m*, defined 
in eq. (0), as a function only of the momentum, i.e., 

. /,n M *( k ) , \ 

<(k) = -j^ 1 - (15) 
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Because of this fact the Hartree terms are exactly calculated, but the Fock ones are 
somewhat approximated as in ref. 0. The exact calculation would imply in considering 
m* as a function of the field operator and this is complicated. However,we still have 
an improvement on the Hartree calculation with the introduction of the exchange terms. 
The nonlinear integral equations for the ZM3 model are the same expressions as obtained 
in the Walecka model with the following modifications: 

1. ) All m s are substituted by m* and all m v are substituted by m*; 

2. ) As a consequence of 1.), the functions Qi(k, q) and $i(k,q) are modified because 
of their dependence on Ai(k,q), which becomes 



where i = o or v . 

For the ZM model just the modifications related to the mass of the sigma meson are 
carried out. All self-energies are evaluated at the self-consistent single-particle energies, 
q° = E(q) and the equations for S s , S° and are solved by a direct iteration procedure. 

We choose to normalize the model parameters using the bulk binding energy and 
saturation density of nuclear matter as usual. As normally done in calculations with the 
Walecka model, we identify the vector meson with the u whose mass is m v = 783 MeV 
and set m s = 550 MeV for the scalar meson mass. For the nucleon mass we take M = 939 
MeV. The energy density for the ZM3 model once the substitution given in eq. (pi) is 
performed yields: 



A(k,q) = k 2 + f + m* 



[E(q) - E{k)f , 



(16) 





dk) + 2g% + Ag 2 v I v , 



(17) 



where the Jj , i = a,v are integrals of the following form 



1 r k F d 3 k r k F d 3 q 



D9(k-q)F i (k,q)H i (k,q) 



(18) 



(2tt) 6 7o E*(k) Jo E*{q) 



with the functions Fi and Hi given by 



F(k, q) = 1/2 - (E(k) - E(q)) 2 D°(k - q), 



(19) 



H a (k,q) = k;q^ + M*(k)M*(q), 
H v (k,q) = k;q*»-2M*(k)M*(q); 



(20) 
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and the D/s are the meson propagators 



D i(k) = T2 72—- • (21) 

To saturate the binding energy per nucleon at —16.15 MeV at the Fermi momentum 
of l.lAfm^ 1 we use g 2 s = 100 and g 2 = 60 for the ZM model. In the ZM3 model, the 
energy saturates at —15.18 MeV at 1.38/m _1 for g 2 = 102 and g 2 = 105. In the Walecka 
model, the energy saturates at —15.75 MeV at 1.31/m _1 for g 2 s = 108.8 and g 2 = 149.4. 

We have then plotted the effective nucleon mass M*/M in function of fcp for the three 
models in fig. 1. One of the reasons for the introduction of the derivative coupling used 
in the ZM and ZM3 models is the small value of the effective nucleon mass at saturation 
density obtained with the Walecka model. We can observe that this problem is corrected 
in the new models. S° and in terms of the momentum also decrease much less in the 
ZM and ZM3 models than in the Walecka model. 

Other important quantities to be analyzed are the scalar and vector potentials and 
the incompressibility. To derive expressions for the potentials, we start from the Dirac 
equation: 

((1 + E v )a ■ k + 7o (M + £ s ) - E°) $ = Eip, (22) 
and after some simple algebraic manipulations, it can be rewritten as 

(* ■ k + 7o(M + 51=-^) + ~ S 1 ° + ^ ) V = (23) 



from where we may define in a natural way 



S = — ; ^ — 24 

1 + T, v K ' 



-S° + EY1" , . 

V = . 25 

1 + E« v ' 



as the scalar potential and 



as the vector potential, where E given by eq . (|14]) can be written as E = Ejp at the 
saturation density. In table 1 these potentials are displayed for the Walecka and the ZM 
models in the mean field and Hartree - Fock calculations. The V — S quantity is related 
to the spin- orbit splitting in finite nuclei and the V + S quantity corresponds to the real 
part of the optical potential for the zero three-momentum. As there is some evidence for 
the strong potentials of the Walecka model, we believe that the ZM3 model can probably 
provide good results for finite nuclei while the ZM potentials are too weak. 
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Concerning the incompressibility K, it has also been calculated for the three models, 
where 

K = 9 Po^-ip=Po (26) 



with po(k) = pBik/kp) 3 - The results are also shown in table 1. According to ref [15 



the expected incompressibility value is K = 210 ± 30 MeV, which means that the results 
obtained for the ZM and ZM3 models are of great improvement in comparison with the 
results coming from the Walecka model. 

To conclude, we would like to comment that we have calculated some relativistic 
features which are important for the understanding of nuclear matter within the context 
of the ZM and ZM3 models. For the present calculation we have taken into account the 
Hartree terms (related to the direct diagrams for the baryon propagator) and Fock - like 
terms (related to the exchange diagrams) with the help of the approximation given by 
eq . fll5j) . We have calculated the effective mass, the incompressibility and the scalar and 
vector potentials and obtained good results out of the ZM3 model. The same treatment 
applied to these models can be extented to finite nuclei calculation in substitution to 
the more complicated relativistic Dirac - Brueckner - Hartree - Fock approach. This 
investigation is under way. 
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Figure Captions 



Fig.l) Effective nucleoli mass as a function of k F . The solid line stands for the Walecka 
model, the dashed line for the ZM model and the dot - dashed line for the ZM3 model. 

Table Caption 

Table 1) Scalar, vector potentials and incompressibility at nuclear matter saturation 
density are shown for the Walecka and the ZM models. MF stands for mean-field and HF 
for Hartree-Fock approximation. 
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models 


S (MeV) 


V (MeV) 


V + S (MeV) 


V - S (MeV) 


K (MeV) 


Walecka - MF 


-431.02 


354.12 


-76.87 


785.18 


550.82 


Walecka - HF 


-458.31 


379.01 


-79.30 


837.32 


585.00 


ZM - MF 


-140.64 


82.50 


-58.13 


223.13 


224.71 


ZM - HF 


-177.5 


109.24 


-68.26 


286.74 


298.47 


ZM3 - MF 


-267.00 


203.71 


-63.28 


470.71 


155.74 


ZM3 - HF 


-264.94 


178.92 


-86.02 


443.86 


161.81 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9603025vl 



